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Abstract 

We have performed full configuration interaction computations of the ground states of 
the molecules Be, BeH2, Li, LiH, B, and BH and verified that the core electrons constitute 
"separated electron pairs." These separated pairs of core electrons have nontrivial structure; 
the core pair does not simply occupy a single spatial orbital. 

Our method of establishing the presence of separated electron pairs is direct and con- 
clusive. We do not fit a separated pair model; we work with the wavefunctions of interest 
directly. To establish that a given group of spin-orbitals contains a quasi-separated pair, 
we verify by direct computation that the quantum state of the electrons that occupy those 
spin-orbitals is nearly a pure 2-electron state. 

Keywords: separated electron pairs, APSG, orbital domains, natural molecular shell 
1. Introduction 

According to separated pair theory , a many-electron wavefunction features an 
ideal "separated electron pair" if and only if some subspace of the 1-electron space, called 
the "orbital domain" of the pair, is occupied by exactly two electrons in a pure state, i.e., 
a state described by a 2-electron wavefunction or "geminal." Where this is approximately 
the case, we will say there is a "quasi-separated electron pair." 
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It has long been known that the electronic structure of the Be atom and LiH 

molecule is very well described by the ansatz of anti-symmetrized products of strongly- 
orthogonal geminals, the APSG ansatz {7, 8, 9]. Though the APSG ansatz does not fit the 
BH molecule very well, the core electrons in BH can still be modeled as a separated pair. 
By 1970, Mehler, Ruedenberg, and Silver 10] had determined that, "The application of 
the separated pair approximation to Li, LiH, B, and BH has shown that the K [i.e., core] 
shell acts as a separated pair and is well described within the context of the separated pair 
approximation." 

The present study confirms this conclusion of Mehler, Ruedenberg, and Silver 10( in 
a direct and compelling way. Ref. IOJ validates the separated pair structure a posteriori, 
by showing how well the separated pair model accounts for correlation energy and predicts 
other observables. In contrast, we discover the structure of the separated electron pairs by 
direct analysis of the wavefunction of interest, not its APSG approximation. 

We first obtain the wavefunction of interest at the full configuration interaction (FCI) 
level of theory. Given a subspace G of the 1-electron space, we derive, from the FCI 
wavefunction, the density matrix Dq on the fermion Fock space over G that describes the 
quantum state of the electrons in G. If Dq is very close to a 2-electron pure state, then G 
may be regarded as the orbital domain of a quasi-separated electron pair. To quantify how 
close Dq is to a 2-electron pure state, we shall use the "2-purity" defined in Section 12.31 

In this way, we have found the orbital domains of the core electron pairs in the ground 
states of Li, LiH, Be, BeH 2 , B, and BH. Each of these core pairs predominantly occupies 
a single spatial orbital core , in the sense that both spin-orbitals CO ret an d <t>corei are 
present in all of the dominant configurations of the FCI wavefunction. The 2-dimensional 
span of the spin-orbitals <p C oret and (f) COT eXi which we call the "trivial core" orbital domain, 
may already be regarded as the orbital domain of a quasi-separated electron pair. However, 
provided that one uses basis sets with enough core orbitals, one observes that the core pair's 
orbital domain involves several orbitals in addition to the trivial core orbitals. To establish 
that those additional orbitals belong in the core orbital domain along with (p core ^ and 
(pcorei, we will show that the state of the electrons in the larger orbital domain, including 
the additional orbitals, is even closer to a 2-electron pure state than is the state of the 
electrons in the trivial core. 
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Finding good orbital domains requires some guesswork. There is a technique for locating 
possible orbital domains which has served us well for the cases treated in the present study, 
though it is not guaranteed to be effective in general. We look for the kinds of statistical 
correlations that an ideal separated pair would entail, as discussed in Section 12.21 Usage 
of this technique is shown by example in Section 13.3.11 but it is not the main point of 
emphasis of the present article. 

The main ideas of this article are that (i) it is feasible to check by direct computation 
whether a low-dimensional subspace of the 1-electron space is the orbital domain of a quasi- 
separated electron pair, and (ii) we have established in this way that the ground states of 
several small molecules actually feature quasi-separated core pairs. 

The rest of this article is organized as follows. Section [2] presents the theoretical back- 



ground of our method, with a short history of separated pair theory in Section 12.11 defini- 



tions and properties of ideal separated electron pairs in Section I2.2[ and a more definitive 
description of our method in Section 12.31 Section |3] demonstrates our method by revealing 
the separated core pairs in ground states of Be, BeH_2, Li, LiH, and B and BH. The Be 
atom is treated in the greatest detail in Section 13.3.11 There is also a brief conclusion 
(Section 0J and an appendix (Appendix A). 



2. Separated electron pairs 

2.1. Brief history of separated pair theory 

The separated electron pair approximation posits the "antisymmetrized product of 
strongly orthogonal geminals" (APSG) ansatz for the wavefunction of the electrons in a 
9|, The APSG ansatz for a 2n-electron wavefunction \1/ is 



molecule 
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* = A 2n (gi ® 92 ® ■ ■ ■ ® 9n) , 



where A2 n denotes an anti-symmetrizing operator and the 2-electron wavefunctions or 
"geminals" gi(x\,X2) are such that 



J ^(1,2)^(1, 2') = 



3 



for i j£ j. The preceding condition of "strong orthogonality" is equivalent 12] to the 
condition that ^ and gj be expressible as linear combinations of Slater determinants of 
spin-orb itals taken from mutually orthogonal subspaces Gi and Gj. 

The APSG ansatz generalizes the Hartree-Fock ansatz, which results when all of the 
geminals in ([1]) are taken to be 2-electron Slater determinants. Thus the APSG ansatz 
variational energy is lower than the Hartree-Fock energy, and at one time it was hoped that 
the APSG ansatz might be broad enough to capture most of the correlation energy l^, 
An early study of the beryllium iso-electronic series found that "uniformly, about 90% of 
the correlation energy is recovered" by the APSG ansatz jsj. Unfortunately, subsequent 
studies found that the separated pair model "works well for LiH but only partially for 
BH" flO|. and demonstrated the "uniform inability of the separated pair wavefunction to 
provide for correlation in the lone-pair geminal in NH" [14J]. It became clear that "even 
an optimum separated-pair function can account for less than half the total correlation 
energy" in larger molecules 

Although APSG was not deemed successful as a simple variational ansatz, much work 
done within the last fifteen years shows that APSG wavefunctions perform very well as 



reference functions for correlated methods 
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Q 3 H, Gi. 



A series of four papers by 



Rassolov, Xu, and Garashchuck has shown that a viable "geminal model chemistry" can 



22|. 



be based on generalized separated pair theory pj, 0, 2l|, 

For more information and historical perspective, we refer the reader to the monograph 
by Hurley [l| that expounds pair theories as a conceptual bridge between Hartree-Fock 
theory and coupled cluster theory, and to a review article by Surjan j^j], which devotes 
several pages to the history of geminal theories. 



2.2. Ideal separated electron pairs 

2.2.1. Definitions for pure and mixed states 

McWeeny 23] generalized the APSG ansatz to antisymmetrized products of strongly 
orthogonal group functions. These are iV-electron wavefunctions of the form 



fk] , 
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where each "group function" fi is an antisymmetric wavefunction of rij variables, Yl n i = N , 
and the group functions are "strongly orthogonal" in the sense that 



J dx 1 f i (l,2,...,n i )f j (l,2',...,n , j ) = 



for % 7^ j. Strong orthogonality of the group functions fx, ■ ■ ■ , fk is equivalent [24( to the 
condition that the group functions g$ can be expanded in Slater determinants taken from 
mutually orthogonal subspaces Gj of the 1-electron spin-orbital space. 

Using this terminology, we can say that a many-electron wavefunction features an "ideal 
separated electron pair" if it equals the antisymmetric product of two strongly orthogonal 
group functions, one of which is a geminal: 

Definition 1. An N -electron wavefunction ^ features an "ideal separated electron pair" 
if there exist strongly orthogonal group functions g(x±, X2) and h(x\, . . . , ccjv-2) such that 

* = A N [g(l,2)h(3,...,N)}. (2) 

If \1/ has the form ([2]), then there exist orthogonal subspaces G and H of the 1-electron 
space such that g is a superposition of 2-particle Slater determinants of spin-orbitals from 
G while h is a superposition of Slater determinants of N — 2 spin-orbitals from H. The 
subspace G belonging to the separated electron pair has been called the pair's "orbital 



domain" [3|, "Arai subspace" [2J, and "carrier space" [25 ]. 

When an ideal separated electron pair with geminal g has orbital domain G, we may 
say that G is occupied by exactly two electrons in the pure state g. This is not just a 
figure of speech. Whatever the iV-electron wavefunction \1/ may be, and whatever the 



subspace of the 1-electron space G may be, the "open subsystem" [26|] consisting of just 
those electrons in G is a many-electron system in its own right, in a quantum state that is 
determined by the state \l/ of the entire molecule. The quantum state of the electrons that 
happen to be in G is typically a mixed state, prescribed formally by a density operator 
on the fermion Fock space over G. In particular, there is typically a "random" number of 
electrons in the subsystem, because electrons are free to move in and out of G. But if G is 
the orbital domain of an ideal separated electron pair, then there are exactly two electrons 
in G, and these two electrons are themselves in a pure 2-electron state described by some 
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geminal. The converse statement holds as well, so that ideal separated electron pairs may 
be characterized as follows: 



Proposition 1. A many-electron wavefunction \1/ features an ideal separated electron pair 
with orbital domain G if and only if the state of the many-electron system consisting of 
just those electrons in G, which is derived from the state \& of the entire system, is a pure 
2 -electron state. 

Definition [1] applies only to "pure" states, those described by many-electron wavefunc- 
tions. Proposition [TJ now allows us to extend the notion of ideal separated electron pairs 
to "mixed" many-electron states, which are described by density matrices: 

Definition 2. A many-electron state features an "ideal separated electron pair" if there 
exists a subspace G of the 1- electron space such that the derived state of the electrons in G 
is a pure 2-electron state. The subspace G is called the "orbital domain" of the separated 
electron pair. 

The formalism wherein the quantum state of the electrons in G is derived from the 
electronic state of the entire molecule is elaborated below in Section 12.31 

2.2.2. Ideal separated electron pairs and natural spin-orbitals 

Eigenvectors of the first-order density matrix (1RDM) derived from a many-electron 
wavefunction are called natural spin-orbitals (NSOs) associated with the wavefunction 27]. 
We shall use the Lowdin normalization of the 1RDM so that its trace equals the number 
of electrons, instead of 1. Then the eigenvalue of the 1RDM corresponding to an NSO is 
called the natural occupation number of 0. We call NSOs of natural occupation number 
"null" NSOs; we are concerned only with non-null NSOs in the following. 

Let g denote a geminal wavefunction. Then a non-null NSO associated with g has even 
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multiplicity as an eigenvector of the 1RDM derived from g. Moreover, there exists 
some orthonormal sequence 0i, 02, • • • of NSOs associated with g such that 

9 = Ci 1 0102 1 + C 2 |0 3 04| + C 3 |0 5 6 | + ••• (3) 

When a many-electron wavefunction \l/ features an ideal separated electron pair with 
geminal g, all NSOs associated with g are also NSOs associated with \1/ with the same 
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natural occupation number. Thus the orbital domain of an ideal separated electron pair is 
spanned by NSOs of the wavef unction. 

The preceding facts imply that the presence of an ideal separated electron pair entails 
certain statistical correlations between occupation events: 

Let \l/ be a many-electron wavefunction and let be a NSO associated with By the 
"multiplicity" of (p we mean its multiplicity as an eigenvector of the 1RDM derived from 
fy. The following four propositions concern NSOs of multiplicity 2. 

Proposition 2. Suppose that features an ideal separated pair with orbital domain G. If 
<p is a non-null NSO of ^ of multiplicity 2 then either (p _L G or (p G G. 

Proposition 3. Suppose that (pi and (p 2 are non-null NSOs of ^> of multiplicity 2, and 
suppose that (pi G G, the orbital domain of an ideal separated pair. Then (pi and (p 2 have 
uncorrected occupations if and only if (p 2 _L G. 

Proposition 4. Suppose that (pi and (p 2 are non-null NSOs of^> of multiplicity 2 that both 
belong to the orbital domain of the same ideal separated pair. If (pi and (p 2 have different 
natural occupation numbers then their occupations are perfectly anti- correlated, in the sense 
that they are never simultaneously occupied by electrons. 

Proposition 5. Suppose that (pi is a non-null NSO of ^ of multiplicity 2 that belongs to 
the orbital domain of an ideal separated pair. If (p 2 is an NSO of ^ that is orthogonal 
to (pi but has the same natural occupation number, then the occupations of (pi and <p 2 are 
perfectly correlated, in the sense that whenever one of the NSOs is occupied then the other 
is also occupied. 

To help us find quasi-separated pairs, we look for correlations between NSOs that 
approximate the ideal behavior described in the propositions. As we demonstrate in Sec- 
tion 13.3. 1\ correlations like those described in the above propositions are easily detected 



when we employ the "marginal-normalized" correlation measure defined in Appendix A 

Correlations like those in the above propositions do not necessarily mean that a sepa- 
rated pair is present. To establish conclusively that there is a separated pair, we proceed 
as described in the next section. 
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2. 3. Quantifying the purity of quasi- separated electron pairs 

A many-electron state features an ideal separated pair if and only if some subspace of 
the 1-electron space is occupied by exactly two electrons in a pure 2-electron state. When 
this situation holds approximately, that is, when some subspace of the one 1-electron space 
is occupied by nearly two electrons in a nearly pure state, we will say that the many-electron 
state features a "quasi-separated" electron pair. 

In the preceding section we noted that the orbital domain of an ideal separated electron 
pair is necessarily spanned by natural orbitals. In contrast, an orbital domain of a quasi- 
separated electron pair need not be spanned by natural orbitals. In Section 13.3.61 we shall 
encounter an example where an orbital domain spanned by canonical Hartree-Fock orbitals 
is better (holds electrons in a state of greater purity) than any orbital domain spanned by 
natural orbitals. 

The principal aim of the present article is to demonstrate the following method of de- 
termining whether a given many-electron state features a quasi-separated electron pair in 
a given orbital domain. Given the many-electron state of interest as a correlated wavefunc- 
tion or mixture of such, and given a subspace G of the 1-electron space, one derives the 
density matrix Dq that describes the quantum state of the electrons in G. One can then 
assess how close Dq is to the density matrix of a pure 2-electron state and decide whether 
or not to regard G as the orbital domain of a quasi-separated electron pair. 

We now explain in more detail how the density matrix Dq is obtained, and how we 
quantify the proximity of Dq to a pure 2-electron state. 

Let G denote a subspace of the 1-electron space H = L 2 (R 3 x {a,f3}). Slater determi- 
nants of n spin-orbitals from G represent configurations of n electrons that all occupy G. 
Let G An denote the span of all such Slater determinants. The fermion Fock space over G is 
used to describe a system of finitely many electrons that all occupy G. This is the Hilbert 
space 

$[G] = C©G©G A2 ©G A3 ©--- (4) 

whose component G An accommodates n-electron states. If dim(G) = d is finite, then 
dim $[G] = 2 d is also finite, because then the d th component of $[G] is 1-dimensional and 
all higher components are O-dimensional. 
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The system consisting of all the electrons in the molecule can be regarded as being 
of two parts: those electrons in G and those in G- 1 , the "orthogonal complement" of G. 



Formally, this is expressed by the isomorphism 29] 



SPeG 1 ] = $[G] <g) ^[G 1 -] . (5) 

The electronic state of a molecule determines the quantum state of its electrons that 
occupy G, as follows. Suppose the electronic state of the entire molecule is represented by 
a density matrix D on the Fock space $[H] = $[G © G ]. Then the density matrix D G 
that represents the derived quantum state of the electrons in G is the partial trace of D 
over the tensor factor ] in (|5]), i.e., 

D G = Tr 5[G x]D . (6) 

The density matrix D G is always block-diagonal with respect to the decomposition 03]) of 
$[G], which we may express by writing 

D G = D G>0 © D G>1 © D G2 © • • • © D Gdim(G) . 

The probability that G contains exactly two electrons equals TrD G 2 . In the event that two 
electrons are found in G, the state of those two electrons is given by the density matrix 
TrD G 2 ^ G ' 2 ' Therefore, there are exactly two electrons in G if and only if TrD G 2 = 1, and 
those electrons are in a pure 2-electron state if and only if the largest eigenvalue of D G 2 
equals 1. 

We denote the largest eigenvalue of D G}2 by \ 2 (G) and call it the "2-purity" of the 
quantum state D G . Thus G is the orbital domain of an ideal separated electron pair if and 
only if X 2 (G) = 1. 

We will use the 2-purity of D G to judge whether or not G may be regarded as the orbital 
domain of a quasi-separated electron pair. Other measures of purity would serve the same 
purpose. The standard measure of purity is the von Neumann entropy. In fact, we originally 
used the von Neumann entropy to measure purity and were led to the same conclusions. 
We prefer to report 2-purity instead of von Neumann entropy because the former suits our 
purpose exactly: the 2-purity of D G indicates whether G contains nearly two electrons in 
a nearly pure state, whereas the von Neumann entropy of D G only quantifies purity and 
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gives no information about particle number. Also, \2(G) is somewhat easier to compute 
than the von Neumann entropy of Dq because it only requires the 2-particle component 

D G ,2- 



3. Results 

3.1. Computational methods 



We used the electronic structure software GAMESS [30[ to obtain the full CI expansions 
of the wavefunctions of interest, using the natural orbitals as reference orbitals. The 
electronic structures of the diatomic molecules LiH, BH, and Li2 were computed with the 
experimental inter-nuclear distances 1.5957A, 1.2324A, and 2.6730A, respectively; and the 
Be-H bond length in the linear molecule H-Be-H was set at 1.34A |3l|. 

The gaussian basis sets we used include the correlation consistent basis sets cc-pVDZ 



321 ] and their "core- valence" extensions cc-pCVDZ |33[. We took the basis set cc-pCVDZ 
for the Be atom from the EMSL Basis Set Exchange 34j, |35|]. We also used some even- 
tempered basis sets of gaussian primitives. 

FCI data obtained from GAMESS became the input for our own programs, which 
we used to derive the density matrices Dq mentioned in Section 12.31 With our current 
algorithms and computing power, computations begin to get quite heavy when dim(G) ~ 
20. 

3.2. Nomenclature for natural orbital domains 

In Section 13.31 we will identify the orbital domains of the separated core pairs in Be, 
BeH 2 , Li, LiH, B, and BH. This section describes the system of nomenclature we use for 
the orbital domains. 

Consider the FCI/cc-pCVDZ ground state wavefunction for the Be atom. We have 
derived and diagonalized its 1RDM. The eigenspaces of the 1RDM are invariant under 
the rotation group, and are in fact irreducibly invariant. We call these eigenspaces of the 
1RDM "irreducible natural molecular shells." We label the irreducible natural molecular 
shells with s,p,d, etc., to signify which irreducible representation of the rotation group 
they support. These eigenspaces are all even dimensional due to spin. The 2-dimensional 
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eigenspaces are numbered ls,2s,3s,4s, etc., in decreasing order of the corresponding eigen- 
values of the 1RDM. Similarly, the 6-dimensional p shells are numbered lp,2p,3p, etc., in 
decreasing order of average occupation. Note that the most occupied natural p shell is 
labelled lp. The numeral 1 in lp is not an energy quantum number; it indexes occupation, 
so lp in this context just means the most occupied irreducible natural p shell. 

In Section l3.3.1l we will establish that the FCI/cc-pCVDZ ground state for the Be atom 
features a quasi-separated pair of core electrons whose orbital domain is not merely the 
natural Is shell, but also contains the irreducible natural 4s and 2p shells. The orbital 
domain is the "direct sum" of the natural Is, 4s, and 2p shells. It will be called the natural 
ls©4s©2p shell. 

The nomenclature for orbital domains is similar for the molecules of our study. For 
example, consider a computed ground state wavefunction for the LiH molecule. We de- 
compose the 1-electron space into irreducible natural molecular shells by diagonalizing the 
1RDM and classifying its eigenspaces by symmetry type a, it, 5, etc. These irreducible 
natural molecular shells are then numbered la,2a,3a, . . .,l7r,27r,37r, . . ., etc., in order of de- 
creasing average occupation number. Having named the irreducible natural shells, larger 
natural molecular shells can be designated using direct sum notation. For example, in 
Section 13.3.41 we will conclude that the FCI/cc-pCVDZ ground state of LiH features a 
quasi-separated electron pair whose orbital domain is the natural lcr©5a"©27r©6<7 shell. 

The ground states of the Li and B atoms are degenerate. In these cases we have 
selected the symmetric mixed ground state for analysis. Since the symmetric ground state 
is invariant under the rotation group, the eigenspaces of its 1RDM can also be classified 
by symmetry as Is, 2s, lp, 2p, etc. 

3.3. Data 

3.3.1. The Be atom 

Quasi-separated electron pairs can be found in FCI ground states of the Be atom 
obtained using various basis sets. Here we report results for the basis sets cc-pVDZ and 
cc-pCVDZ, as well as an even-tempered basis set. 

First we present the results for the FCI/cc-pCVDZ ground state. Some of the data for 
this wavefunction is presented in Table [TJ The natural occupation numbers of the NSOs 
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are listed in the Occupation column. Note that these are the average occupations of the 
natural spin-orbitals, not the natural shells. For example, to get the average occupation 
of the irreducible natural lp shell, the spin-orbital occupation 0.02979 in Tabled] must be 
multiplied by 6, since the lp shell is spanned by six lp NSOs. 

The main part of Table [1] presents the values of the marginal-normalized correlation 
between NSO occupations. Marginal-normalized correlation, denoted by k, is defined in 



Appendix A It is designed to detect the correlations expressed in Propositions [3J15] of 



Section 12.21 For uncorrelated occupations, as in Proposition [3] the value of k is 0. For 
perfectly anti-correlated occupations, as in Proposition HJ the value of k is — 1. For perfectly 
correlated occupations, as in Proposition El the value of k is 1. Values of n are given to 
one decimal place. An entry -0 in the table indicates that k is negative, and an entry 1 or 
-1 in bold face indicates that 0.99 < \k\ < 1. 

The values of k given in Table [1] reveal the quasi-separated electron pair structure. The 
2s, lp, 3s, and Id NSO occupations appear to be very highly correlated with one another 
but not with the Is, 4s, and 2p NSO occupations. We therefore suppose - and will soon 
confirm - that the 2s, lp, 3s, and Id NSOs span the orbital domain of a quasi-separated 
electron pair, for this hypothesis is consistent with Propositions [2] - [5j The Is, 4s, and 2p 
NSO occupations also appear to be highly correlated with one another, consistent with 
the hypothesis that they span the orbital domain of a quasi-separated electron pair. The 
3p NSOs appear to be moderately correlated with NSOs from both presumptive separated 
pairs, but they are very lightly occupied and may be neglected. 

We now argue that the core electrons form a quasi-separated pair whose orbital domain 
is the natural ls©4s©2p shell. Of course, the natural Is shell by itself may be also regarded 
as the orbital domain of a quasi-separated pair, for it is occupied by almost 2 electrons in 
the pure Slater determinant state |ls Q ls/3|. However, including the 2p NSOs in the orbital 
domain of the core pair increases the purity of the state, and including the 4s NSOs as well 
increases the purity even further. We can quantify this increase in purity by measuring 
the "2-purity" A 2 defined in Section 12.31 

Recall that \2{G) denotes the largest eigenvalue of the 2-particle component of the 
density matrix Dq that represents the state of the electrons in G, so that A2(G) = 1 if and 
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0.00007 


-0 


-1 


-0.8 


-0.8 


-1 


-0.9 


-0.1 











-1 


-1 


-1 


-1 


1 


-0.3 








3p 


0.00001 


-0.4 


-0.4 


-0.7 


-0.8 


0.3 











-0.2 


0.2 











-0 


-0.1 


0.3 


-0.9 


-0.9 


3p 


0.00001 


-0.4 


-0.4 


0.3 


-0.9 


-0.7 











0.1 























0.3 


-0.9 


3p 


0.00001 


-0.4 


-0.4 


-0.9 


0.3 


-0.8 








0.1 





























0.3 



Table 1: Marginal- normalized correlations between NSO occupations of the FCI/cc-pCVDZ ground state 
of the Be atom. Correlations between NSOs of opposite spin appear on and below the main diagonal; 
correlations between NSOs of the same spin appear above the main diagonal. 
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only if G is the orbital domain of an ideal separated electron pair. We observe that 

A 2 (ls) = 0.9985 < A 2 (ls © 4s) = 0.9987 < A 2 (ls © 2p) = 0.9992 

< A 2 (ls © 4s © 2p) = 0.9994 

and conclude that the natural ls©4s©2p shell is the orbital domain of a quasi-separated 
pair of core electrons. 

Next we present the results for the smaller basis set cc-pVDZ. The marginal-normalized 
correlations tabulated in Table [2] suggest that the 2s©lp©3sffild shell may be the orbital 
domain of a quasi-separated pair. Indeed it is, for its 2-purity is 0.9999. The complementary 
shell, the natural ls©2p shell, has the same 2-purity a priori and may also be regarded as 
the orbital domain of a quasi-separated pair. 

However, this core pair is virtually trivial: it does not amount to much more than a 
doubly occupied spatial Is orbital. The natural ls©2p shell has 2-purity 0.9999, but the 
trivial core orbital domain (the natural Is shell by itself) actually has slightly higher 2- 
purity. Moreover, the data in Table [2] does not bear the signature of a distinctly separated 
core pair. According to Proposition [5], if the Is NSOs belonged to the orbital domain of an 
ideal separated pair of core electrons, then the marginal-normalized correlation k between 
the two Is NSOs would equal 1, but the actual value is only 0.4. From Proposition [3j 
we would also expect the Is NSO and 2s NSO occupations to be nearly uncorrelated, but 
k = 0.6 between Is and 2s NSOs of opposite spin. The preceding considerations suggest 
that, with cc-pVDZ, the core electron pair is not distinctly separated. 

Apparently, the structure of the quasi-separated core pair is better revealed when we 
use the larger "core- valence" basis set cc-pCVDZ. We believe the reason is that the larger 
basis set includes more core orbitals and is therefore better able to reveal a core pair that 
is actually present in the "exact" wavefunction, i.e., in the infinite basis set limit. To 
test this, we have tried using large even-tempered basis sets. These trials have repeatedly 
confirmed that the natural ls©4s©2p shell deserves to be regarded as the orbital domain 
of a quasi-separated core pair. 

For example, here are the results for an even-tempered basis set with 10 primitive 
gaussian s-orbitals and 4 primitive gaussian p-orbitals. We took the exponents of the 



primitives from Table I of 



36j. The exponents of the s-orbitals were those of the N s = 10 
14 



Occupation 


Is 


2s 


IP 


lp 


IP 


3s Id Id Id 


Id 


Id 


2p 


2p 


2p 


Is 


0.99996 


0.4 


0.1 


-0.1 


-0.1 


-0.1 


-0.1 


-0 


-0 


-0 


-0 


-0 


-0.6 


-0.6 


-0.6 


2s 


0.90798 


0.6 


1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-0.4 


-0.4 


-0.4 


lp 


0.03010 


-0.1 


-1 


1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-0.6 


-0.1 





lp 


0.03010 


-0.1 


-1 


-1 


1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-0.6 





-0.1 


lp 


0.03010 


-0.1 


-1 


-1 


-1 


1 


-1 


-1 


-1 


-1 


-1 


-1 





-0.6 


-0.6 


3s 


0.00138 


-0.1 


-1 


-1 


-1 


-1 


1 


-1 


-1 


-1 


-1 


-1 











Id 


0.00007 


-0 


-1 


-1 


-1 


-0.9 


-1 


1 


-1 


-1 


-1 


-1 











Id 


0.00007 


-0 


-1 


-1 


-0.9 


-0.9 


-1 


-1 


1 


-1 


-1 


-1 











Id 


0.00007 


-0 


-1 


-0.9 


-1 


-0.9 


-1 


-1 


-1 


1 


-1 


-1 











Id 


0.00007 


-0 


-1 


-0.9 


-0.9 


-1 


-1 


-1 


-1 


-1 


1 


-1 


-0.3 








Id 


0.00007 


-0 


-1 


-0.9 


-0.9 


-1 


-1 


-1 


-1 


-1 


-1 


1 


-0.6 








2p 


0.00001 


-0.4 


-0.6 


-0.7 


-0.7 


0.5 














-0.1 


-0.5 


0.5 


-0.9 


-0.9 


2p 


0.00001 


-0.4 


-0.6 


0.2 


0.3 


-0.7 























0.5 


-0.9 


2p 


0.00001 


-0.4 


-0.6 


0.3 


0.2 


-0.7 


























0.5 



Table 2: Marginal-normalized correlations between NSO occupations of the FCI/cc-pVDZ ground state 
of the Be atom. Correlations between NSOs of opposite spin appear on and below the main diagonal; 
correlations between NSOs of the same spin appear above the main diagonal. 
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Occupation 


Is 


2s 


IP 


lp 


IP 


3s 


4s 


2p 


2p 


2p 


3p 


3p 


3p 


5s 


4p 


4p 


4p 


Is 


0.99830 


0.9 





-0 


-0 


-0 


-0.1 


-0.8 


-1 


-1 


-1 


-0.5 


-0.5 


-0.5 


-0.6 


-0.3 


-0.3 


-0.3 


2s 


0.90742 


0.1 


1 


-1 


-1 


-1 


-0.9 


-0.2 


-0.1 


-0.1 


-0.1 


-0.7 


-0.7 


-0.7 


-0.5 


-0.8 


-0.8 


-0.8 


lp 


0.03025 


-0 


-1 


1 


-1 


-1 


-0.9 


-0.2 


-0.1 


-0 





-0.6 


-0.9 


-0.9 


-0.4 





-0.9 


-0.9 


lp 


0.03025 


-0 


-1 


-1 


1 


-1 


-0.9 


-0.2 


-0.1 





-0 


-0.9 


-0.9 


-0.6 


-0.4 


-0.9 


-0.9 





lp 


0.03025 


-0 


-1 


-1 


-1 


1 


-0.9 


-0.2 





-0.1 


-0.1 


-0.9 


-0.6 


-0.9 


-0.4 


-0.9 





-0.9 


3s 


0.00192 


-0.1 


-0.9 


-0.8 


-0.8 


-0.8 


0.9 























-0.3 











4s 


0.00050 


-0.8 


-0.2 


-0.2 


-0.2 


-0.2 


0.4 


0.6 























-0.4 


-0.4 


-0.4 


2p 


0.00032 


-0.9 


-0.1 


-0 


-0.1 


0.1 








0.8 


-1 


-1 





0.2 








-0.3 





-0.4 


2p 


0.00032 


-0.9 


-0.1 





0.1 











-1 


0.8 


-1 








0.2 





-0.4 


-0.3 





2p 


0.00032 


-0.9 


-0.1 


0.1 





-0.1 








-1 


-1 


0.8 


0.2 














-0.4 


-0.3 


3p 


0.00004 


-0.5 


-0.2 


0.2 


-0.9 


-0.9 














0.3 


0.2 


-0.9 


-0.9 





0.1 








3p 


0.00004 


-0.5 


-0.2 


-0.9 


-0.9 


0.2 








0.3 








-0.8 


0.2 


-0.9 








0.1 


(j 


3p 


0.00004 


-0.5 


-0.2 


-0.9 


0.2 


-0.9 











0.3 





-0.8 


-0.8 


0.2 











0.1 


5s 


0.00002 


-0.5 


-0.5 


-0.2 


-0.2 


-0.2 


0.1 


0.2 




















0.6 











4p 


0.00001 


-0.1 


-0.8 





-0.9 


-0.9 














0.1 


0.2 











0.5 


-0.9 


-0.9 


4p 


0.00001 


-0.1 


-0.8 


-0.9 


-0.9 











0.1 











0.2 








-0.9 


0.5 


-0.9 


4 P 


0.00001 


-0.1 


-0.8 


-0.9 





-0.9 











0.1 











0.2 





-0.9 


-0.9 


0.5 



Table 3: Marginal-normalized correlations for the FCI ground state of the Be atom with an even- 
tempered basis set. Correlations between NSOs of opposite spin appear on and below the main diagonal; 
correlations between NSOs of the same spin appear above the main diagonal. Ten NSOs of occupation 
less than 0.000005 are omitted. 

series given for the Be atom, and the p-orbitals' exponents were those of the N p = 4 series 
given for the B atom, scaled by || = |^ to account for the main effect of the two atoms' 
differing nuclear charges. Table [3] shows the marginal-normalized correlations. Again the 
correlations between the Is and 2p NSO occupations stand out. Again, the 2-purity of the 
natural ls©4s©2p shell is greater than that of the natural Is shell alone: 

A 2 (ls) = 0.9981 < A 2 (ls © 4s) = 0.9982 < A 2 (ls © 2p) = 0.9989 

< A 2 (ls © 4s © 2p) = 0.9990 . 

We tried several other even-tempered basis sets and got similar results. 

Our analysis shows that the Be atom ground state features a quasi-separated core 
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pair in the natural ls©4s©2p shell. The separated core pair we observe appears to have 
the structure predicted 



numbers in Table XI of 



by the separated pairs analysis of Q, judging by the occupation 
or Table 2.6 in 



The character of the core pair has been evident since at least 1960, when Linderberg 



and Shull H| reported that the inner pair of electrons in a 3- or 4-electron atom is very 
similar to the corresponding electron pair in the He iso-electronic series, while they judged 
correlation to be negligible between the shell of the core and the shell of the valence in Be. 
They wrote, "We may therefore conclude that the inner pair of electrons in these atoms 
[Li and Be] is very similar to the pair in the He-like series. This is a not unexpected result 
since a spherical shell of charge contributes nothing to the potential within the sphere, and 
the outer electrons in Li and Be-like atoms are largely in regions of space clearly outside 
the domain of the inner shell." 

We have also investigated ground states of B + , C 2+ , and N 3+ , which are iso-electronic 
with Be. When we use cc-pCVDZ or large even-tempered basis sets to compute them, 
the ground states of these ions are generally found to possess quasi-separated core pairs in 
their natural ls©2p shells. By comparing 2-purities, we can clearly see that the orbitals 
homologous to the natural 4s orbitals in the Be core do not belong in the C 2+ andd N 3+ 
core orbital domains. In B + , however, there is a case to be made for inclusion of the 4s 
NSOs in the core domain. Using cc-pVDZ, the 2-purities do increase slightly when the 4s 
NSOs are included: 

A 2 (ls) = 0.99903 < A 2 (ls © 4s) = 0.99906 < A 2 (ls © 2p) = 0.99947 

< A 2 (ls © 4s © 2p) = 0.99949 . 

On the other hand, a slightly negative effect of including the 4s orbitals is observed when 
we use the von Neumann entropy S(Dq) = — Tr(ZAs* log 2 Dq) to measure the purity of 
Dq- Including the 4s NSOs in the orbital domain G increases the entropy (and thereby 
decreases the purity) of D G : 

S(D ls ) = 0.0118 < S(D ls(B43 ) = 0.0129 < S(D ls(B2p ) = 0.0078 

< S( J D 1 se4se2 P ) = 0.0083 . 

Thus, the 4s NSOs appear to evaporate out of the core pair's orbital domain as we 
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move up the Be-isoelectronic series from Be to N 3+ . An anonymous referee remarks that 



this phenomenon probably reflects the fact [37|, |38|, |39_|, I40J that the ratio of "radial correla- 



tion" to "angular correlation" 4l| decreases along the He iso-electronic series, for the core 
electron pair in a Be-isoelectronic atom is expected to behave much like the two electrons 
in the corresponding He-isoelectronic atom. 

3.3.2. BeH 2 

The FCI/cc-pCVDZ ground state of BeH2 features a quasi-separated pair of core elec- 
trons in its natural I<7 5 ffi5cr g ©4cr u ffi37r n shell: 

\ 2 (la g ) = 0.9984 < A 2 (la 9 © 5a g © Aa u © 3tt u ) = 0.9992. 

The Icr g ffi5(7 5 ©4c7 n ffi37r u shell of the core electrons in BeH 2 corresponds closely to the 
ls©4s©2p shell of the core electrons in the Be atom, with the 2p shell of Be splitting into 
the 4cr u and 3ir u shells of BeH 2 . 

3.3.3. The Li atom 

The ground state of the Li atom is degenerate 2 S. The ground state subspace is spanned 
by two wavefunctions ty a and ^ p of opposite spin. We consider the mixed ground state 
ll^aX^oJ + II^X^I and label the eigenspaces of its 1RDM as described in Section I3T21 

In the FCI/cc-pCVDZ ground state, 

A 2 (ls) = 0.9971 < A 2 (ls © 3s) = 0.9983 < A 2 (ls © lp) = 0.9985 

< A 2 (ls © 3s © lp) = 0.9997. 

With the smaller basis set cc-pVDZ, the trivial core has the greatest purity, with A 2 (ls) = 
0.9999. However, when the basis set cc-pVDZ is de-contracted, and we make a basis set 
out of all of its primitives, the natural ls©3sffilp shell reemerges as the orbital domain of 
the core. 

3.3.4. LiH 

The FCI/cc-pCVDZ ground state of LiH features a quasi-separated core pair located 
in the natural lcr©5<x©27r©6cr shell: 

A 2 (la) = 0.9971 < A 2 (la © 5a © 2tt © 6a) = 0.9994. 
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This orbital domain is closely related to the orbital domain of the core electrons in the 
Li atom, the natural ls©3s©lp shell: the natural lp shell in Li splits into the 2n and 6a 
shells in LiH. 

We do not detect an interesting quasi-separated core pair when we use the smaller 
basis set cc-pVDZ. The orbital domain of the quasi-separated core pair in the FCI/cc- 
pVDZ ground state of LiH amounts to little more than the natural la shell, which already 
has 2-purity 0.9999. Including other orbitals brings the average occupation closer to 2 
electrons but does not increase purity. However, when we use large even-tempered basis 
sets, we again discover a quasi-separated core pair with orbital domain I<7©5crffi27rffi6cr. 

3.3.5. The B atom 

The ground state of the B atom is 2 P, so the ground state subspace is 6-dimensional, 
spanned by six orthonormal wavefunctions ^ z a , We choose a definite 

ground state to work with. Any wavefunction in span{\l/^, . . . , represents a pure 
ground state, and any density matrix that is spectrally supported on span{\l/^, . . . 
represents a mixed ground state. Our choice is the symmetric mixed state 

i\k)(k\ + + i\k)(k\ + i\nm\ + §i«?x*si + hnm\ • w 

This state is invariant under the group action, so we may name the irreducible natural 
shells just as described in Section 13.21 

In the FCI/cc-pCVDZ ground state of the form ([7]), 

A 2 (ls) = 0.9990 < A 2 (ls © 4s) = 0.9992 < A 2 (ls © 3p) = 0.9994 

< A 2 (ls ©3p© 4s) = 0.9996 . 

We conclude that the symmetric mixed state (|7|) features a quasi-separated core pair in its 
natural ls©3p©4s shell. 

3.3.6. BH 

The FCI/cc-pCVDZ ground state of BH features a quasi-separated pair of core electrons 
in its natural lcr©47r©10cr shell: 

A 2 (1ct) = 0.9990 < A 2 (la©47r) = 0.9991 < A 2 (la © 4tt © 10a) = 0.9992. 
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As mentioned in Section 12. 3[ the "best" orbital domain of a quasi-separated electron 
pair need not be a natural molecular shell. For the FCI/cc-pCVDZ ground state of BH, 
the best orbital domain we have found is a Hartree-Fock (HF) molecular shell, spanned by 
canonical HF orbitals. This orbital domain K is the direct sum of three irreducible HF a 
shells and one irreducible HF n shell. The 2-purity of the state Dk is 0.9996. 

The orbital domain K would seem to be a slight perturbation of the orbital domain 
found for the core pair of the B atom in the preceding section. We have not computed the 
overlaps, but occupation numbers tell the story. The irreducible a shells in K have average 
occupations 1.99812, 0.00067, and 0.00032, and the 7r shell has average occupation 0.00065. 
The a shells with average occupation 1.99812 and 0.00067 correspond to the natural Is and 
4s shells of the symmetric mixed state (ED) of the B atom, which have average occupations 
1.99817 and 0.00065, respectively. The irreducible natural 3p shell in the orbital domain 
of the B atom core seems to have split into the n shell of average occupation 0.00065 and 
the a shell of average occupation 0.00032, for 0.00065 + 0.00032 = 0.00097 is very close to 
the average occupation of the B atom's natural 3 p sh ell, i.e., 0.00098. This is all consistent 
with the results of the separated pair analysis of 10( concerning the "transferability of the 
K [core] shell pair in B and BH." 

We do not find a nontrivial core orbital domain when we use the smaller basis set cc- 
pVDZ. However, when this basis set is de-contracted, and we use the basis set consisting 
of all the primitive gaussian orbitals that contribute to cc-pVDZ, we do find that 

X 2 (la) = 0.9987 < A 2 (l<r © 4vr © 10a) = 0.9988 < A 2 (lcr © 4tt) = 0.9989 . 
4. Conclusion 

The presence of quasi-separated electron pairs can be established conclusively by direct 
computation when a sufficiently correlated representation of the many-electron state of 
interest is available. It is convenient to work with the wavefunction in CI form, but it need 
not be the full CI wavefunction. We limited the present study to FCI wavefunctions in 
order that our specific conclusions concerning the ground states of Li, LiH, Be, BeH2, B, 
and BH would be as compelling as possible. In future work, we may apply our method to 
larger molecules by working with coupled cluster wavefunctions. 
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Using our method, we have shown that FCI ground states of Li, LiH, Be, BeH-2, B, 
and BH all appear to feature interesting quasi-separated pairs of core electrons, provided 
large enough basis sets are used. In each case, we have identified the main components of 
the core pair's orbital domain. Our results are consistent with the predictions of separated 
pair theory jl, 6, 10 1. 
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Appendix A. Marginal-normalized correlation 

In this section we define the "marginal- normalized correlation" k that we used in Sec- 
tion [3]3J] to test for statistical correlations like those in Propositions [3]- [5] of Section l2~2l 

Suppose that 0i and 02 are two NSOs associated to a many-electron state, and let 
G = span{0!, cf>2 }• The quantum state Dq of the electrons in G is characterized by three 
numbers, namely, the probability n\ that 0i is occupied, the probability n 2 that 02 is 
occupied, and the probability p that 0i and 02 are occupied simultaneously. The probability 
that (pi and 02 are both unoccupied is then 1 — ri\ — ri2 + p, and it can be shown that 
there are no "off-diagonal" correlations. The marginal- normalized correlation between two 
NSOs is simply a function of ni, n 2 , and p. 

Let us examine the statistical correlation between the occupations of two NSOs from our 
data set. We have selected two NSOs whose occupation probabilities are 0.9288 and 0.0010, 
rounded to 4 decimal places. The probability that both NSOs are occupied simultaneously 
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is 0.0003. The joint occupation probabilities can be arranged in a 2 x 2 table: 





vacant 


occupied 




vacant 


0.0705 


0.0007 


0.0712 


occupied 


0.9285 


0.0003 


0.9288 




0.9990 


0.0010 





The numbers on the right and bottom are the marginal sums, which are the occupation 
and vacancy probabilities for the two NSOs. 

From the preceding probabilities, we can compute measures of correlation. We may 
assign numerical values to "vacancy" and "occupation" and compute the familiar (Pearson) 
correlation coefficient of the resulting random variables (it does not matter which numerical 
values one chooses; the correlation coefficient will come out the same because there are 
only two events: occupation and vacancy). For the probabilities in table (1 A. 1 [) above, 
p « -0.077. 

The correlation coefficient p does not suit our purposes because it does not take the 
marginal probabilities into account. For example, a table of joint probabilities with the 
same marginals as table ( lA.ljl above can only have p between —0.114 and 0.009 (rounding 
to 3 places). These extremes are realized by the tables 



0.0702 


0.0010 


0.0712 




0.0712 


0.0000 


0.0712 


0.9288 


0.0000 


0.9288 


and 


0.9278 


0.0010 


0.9288 


0.9990 


0.0010 






0.9990 


0.0010 





We prefer to use a measure of correlation that depends on the marginals, but which can 
attain a minimum of —1 and a maximum of +1. The variable p in a table of probabilities 



1 — ni — n 2 + p 


n 2 - p 


1 — ni 


ri\ — p 


p 


m 


1 — n 2 


n 2 





can range between max{0,n! + n 2 — 1} and min{ni, n 2 }, and independence occurs when 
p = n\U 2 . The greatest that \p — n\n 2 \ can be is 

min{n!, n 2 } — ri\n 2 if p > n\n 2 

n\n 2 — max{0, n\ + n 2 — 1} if p < n\n 2 ; 
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accordingly, we define the marginal-normalized correlation k (m, n 2 , p) to be 

p - n x n 2 



min{ni, n 2 } — n\n 2 

and 



when p > n\n 2 



P ~ n x n 2 

when p < n\n 2 . 



n\n 2 — max{0, n\ + n 2 — 1} 
The probability tables in (1A.2|) thus have k = —1 and +1. The occupations of the two spin 
orbitals that are tabulated in flA.lj) have a small correlation coefficient p ~ —0.08, but their 



marginal- normalized correlation k « —0.68. This value of k indicates that, for two spin- 
orbitals whose occupation probabilities happen to be 0.9288 and 0.0010, the occupations of 
these two spin orbitals are markedly anti-correlated. 

The value of k is very sensitive to the precision of the occupation probabilities. Had we 
rounded the occupation probabilities to 5 decimal places instead of 4, k would have come 
out to be —0.66 instead of —0.68. Fortunately, we do not require too much precision in k, 
as we usually round it to one decimal place. 
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